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Precise measurement of the refractive index of stable silicon isotopes 28 Si, 29 Si, 30 Si sin gle crystals 
with enrichments above 99.9 at.% and a silicon single crystal nat Si of natural isotopic composition 
is performed with the Fourier-transform interference refractometry method from 1.06 to more than 
80 um with 0.1 cm - resolution and accuracy of 2 x 10" 5 ... 1 x 10" 4 . The oxygen and car- 
bon concentrations in all crystals are within 5 x 10 15 cm -3 and the content of metal impurities is 
10 -5 . . . 10 -6 at.%. The peculiar changes of the refractive index in the phonon absorption region 
of all silicon single crystals are shown. The coefficients of generalized Cauchy dispersion function 
approximating the experimental refractive index values all over the measuring range are given. The 
transmission and Raman spectra are also studied. 
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I. INTRODUCTION 

Production and study of the properties of isotopically 
enriched substances are one of the crucial problems of 
contemporary material science involving in itself the in- 
creasing number of researchers. Owing to the achieved 
level of purity of many substances, their properties are 
substantially determined by the isotopic composition. 

Development of the technology of high-purity monoiso- 
topic single crystals production allows one to study the 
influence of the isotopic composition on physical prop- 
erties, in particular on the electronic, phonon, and Ra- 
man spectra and the refractive index dispersion. This 
is important both for verification of existing theoretical 
models and filling them with more correct values of used 
parameters, and for calculation of parameters of next- 
generation devices developed on their basis. 

Crystalline silicon is known to be the basic material of 
microelectronics and solar energetics. Besides, it is also 
an attractive optical material with high transmission in 
the near- to mid-IR spectral range. In recent years sil- 
icon photonics grows rapidly, which is expected to find 
application in future optical fiber communication and in- 
formation transfer systems. Monoisotopic silicon is con- 
sidered to be the basic component of quantum computers 

Formerly silicon refractive index was measured repeat- 
edly in various spectrum ranges by many laboratories 
and companies (see e.g. [i~7l - l2ll p. The papers [l7], HI] 
seem to be the first in this field. The analysis of more 
than 50 papers made in 1980 [19| showed the data avail- 
able to be scattered widely, but allowed to derive certain 
averaged spectral dependence for silicon refractive index 
which is rather often referred to. Unfortunately, in the 
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majority of published papers devoted to silicon refractive 
index measuring, the crystals used for this purpose are 
characterized neither in purity, nor in conductivity, both 
influencing the measurement results. 

In this paper we present the results of measuring 
the room-temperature transmission spectra in the short- 
wavelength edge and phonon absorption regions, the Ra- 
man spectra excited by the 514.5 nm Ar laser line and 
the refractive index spectral dependence in the near-, 
mid-, and far-IR spectral ranges in silicon single crystals 
with natural isotopic composition, nat Si, and silicon iso- 
topes 28 Si, 29 Si, 30 Si with the basic isotope content above 
99.9 at.% prepared by the same technology. A short note 
on fabrication of these single crystals and refractive in- 
dex measurements in the 1.06 ... 25.5 um range has been 
published recently Q. 



II. PREPARATION OF SINGLE CRYSTALS 
AND SAMPLES FOR MEASUREMENTS 



Polycrystalline silicon was prepared by thermal decom- 
position of monosilane enriched in one of the silicon iso- 
topes as described in [H, H| . The decomposition was pro- 
ceed on a graphite substrate surface at 800 °C. After a 
few mm thick layer was deposited, polycrystalline silicon 
was separated from the graphite substrate, remelted to 
form a rod, and purified by 10 float zoning passes. Then 
this monoisotopic silicon rod was used as a substrate to 
deposit the polycrystalline silicon layer till an ingot of a 
required diameter was obtained. 

To minimize the effect of isotopic dilution by a seed- 
ing material, a special growth procedure [j| was used. 
Preliminarily a single-crystal seed with the isotopic com- 
position gradient in length was grown by drawing at the 
end of mono-isotopic silicon ingot of a cylindrical site 
with a diameter not exceeding that of a seed from sili- 
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Table I: Results of mass spectrometry determination of 
isotopic composition of produced Si single crystals 



Single 
crystal 


Isotopic composition, at.% 


28 Si 


29 Si 


30 Si 


28 Si 


99.9934 


0.00637 


0.00023 


29 Si 


0.026 


99.919 


0.055 


30 Si 


0.005 


0.021 


99.974 



con of natural isotopic composition. The composition of 
grown single crystal seed was almost identical to the iso- 
topic composition of polycrystalline ingot which the sin- 
gle crystal was grown from. Single crystals were grown in 
the [100] direction by float zoning in a high-purity argon 
atmosphere. The same procedure was used to prepare a 
single crystal of natural isotopic composition, nat Si: 28 Si 
(92.23%), 29 Si (4.68%), and 30 Si (3.09%) @j. 

The isotopic composition of the crystals was deter- 
mined by mass spectrometry Q with an accuracy bet- 
ter than 0.01 at.% (Table HJ. Comparison with the data 
reported in 0, Q indicates that our crystals are more iso- 
topically enriched and offer higher chemical purity and 
resistivity. 

Structural perfection of the crystals was evaluated us- 
ing X-ray diffraction and selective etching. No orienta- 
tion disorder, low- angle boundaries or twinning regions 
were detected in the single crystals grown. The half- 
width of their rocking curves was 11 . . . 12". The room- 
temperature resistivity of the crystals was found to be 
100-200 x cm. According to mass spectrometric mea- 
surements, the contents of 72 impurities was below the 
mass spectrometry detection limit (< 10~ 5 . . . 10 -6 at.%) 

2 ... 3 mm-thick plates were cut perpendicularly to the 
growth axis of the cylindrical samples prepared. Then 
the plates were ground and polished with various flatness 
degrees for measuring the transmission spectra without 
interference at 2 cm" 1 resolution (to analyze the presence 
of impurity absorption bands and intrinsic bands shift in 
the two-phonon absorption range) and with interference 
at 0.1 cm -1 resolution (to measure the refractive index 
spectral dependence). Samples for the refractive index 
measurements had the surface flatness better than 1/10 
of Newton's fringe (A = 546 nm), parallelism of sides 
better than 1", and the surface finish class of 60/40. 



III. TRANSMISSION SPECTRA 

The transmission spectra of crystalline silicon without 
interference were measured on Lambda 900 spectropho- 
tometer in a short- wavelength range (0.9 . . . 1.6 pm) and 
on Bruker IFS-113v Fourier-transform vacuum spectrom- 
eter in a longer wavelength range (1...90 pm). The 
analysis on the presence of impurity absorption bands 
in the transmission spectra measured in conformity with 



the ASTM standards F-1188 and F-1391 [III proved oxy- 
gen and carbon contents to be less than 5 x 10 15 cm~ 3 
in all the samples under investigation. 

In the 1 to 3 mm thickness crystals the tendency to 
the transmission edge shift to higher frequencies with iso- 
tope atomic mass increasing was observed in the short- 
wavelength transmission edge region at room tempera- 
ture. However, the shift value of < 1 meV/a.m.u. was 
comparable with the experimental error. 

The M -1 / 2 dependence of the band gap, E g , was 
shown in [l2| to be more probable for T < 50 K, since 
due to lattice expansion and electron-phonon interac- 
tion Eg decreases linearly with T growing in the room- 
temperature region (~ 294 K). Measurement at 16 K of 
the transmission spectra of our silicon samples showed 
that with atomic mass increasing the absorption edge 
shifts to the high-energy region. The shift was found to 
be T (1-8 ±0.3) meV in 28 Si and 30 Si relative to 29 Si. 
The indirect transition energy determined from the opti- 
cal absorption edge was found to depend on the average 
atomic mass as Af -1 / 2 . The shift ~ 1.5 meV measured 
in the present work at 16 K turned out to be somewhat 
higher compared with 1 meV derived from the lumines- 
cence spectra in 

In the phonon absorption region (5 ... 30 pm) a shift of 
all absorption bands to longer wavelengths was observed, 
which had been seen previously in isotopically enriched 
silicon single- crystal samples with a smaller isotopic en- 
richment 11311. 



IV. RAMAN SPECTRA 

The Raman spectra of all produced single crystals were 
measured on a T64000 triple spectrograph in a backscat- 
tering geometry at room temperature. The 514.5-nm ar- 
gon laser line was used as an incident radiation. The scat- 
tered radiation was detected by a liquid nitrogen-cooled 
CCD matrix. The spectra were gathered in 0.7 cm -1 
steps. Spectrograph was calibrated using radiation in Ar 
laser and Hg lamp lines. 

The Raman band corresponding to first-order scatter- 
ing by Ti5 symmetry phonons was observed in spectra of 
all the single crystals (showed as circles in Fig. Ha]) . The 
band turned out to be well fitted by a Lorentzian with 
width 3.3 .. . 3.4 cm" 1 for 28 Si, 29 Si, 30 Si and 3.7 cm- 1 for 
nat Si (solid lines in Fig. Hal) . This is an evidence of high 
structural perfection of the mono-isotopic single- crystal 
samples. 

The peak position, f2, of this Raman band is plotted 
vs. atomic mass, Mi, in Fig.[Tb]for all silicon single crys- 
tals. The dashed line represents the best fit of the exper- 
imental phonon frequency values for all the single-crystal 
samples by the VL = 2731.9 {Mi [a.m.u.) 1 ^ 2 cm -1 depen- 
dence valid in harmonic approximation. The root-mean- 
square deviation, SD, of experimental f2 points from the 
fitted line turned out to be small for the mono-isotopic 
samples. Again this is an indirect evidence of high struc- 
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(a) 



(b) 



Figure 1: (a) Fundamental Raman band of 28 Si, 29 Si, 30 Si, and nat Si (natural isotopic composition) single crystals. 
|(b)| Peak position f2 of the Raman band vs. isotopic crystal composition 



tural perfection of those. The Raman band peak position 
for nat Si sample falls slightly above the fitted line, the 
band being wider than in mono-isotopic samples. This 
is attributable to both anharmonicity related to different 
atomic masses of isotopes and to worse structural perfec- 
tion of nat Si crystal in comparison with the mono-isotopic 
ones (see further). 

V. REFRACTIVE INDEX 

A. Measurements 

To measure refractive index we used the improved 
method of interference refractometry fT^ - fl6| . In this 
method firstly the order, m, of the interference maximum 
is determined with the absolute accuracy from the peak 
position, v m , in the interference transmission spectra of 
plane-parallel plates cut from the studied material, and 
then the refractive index value, n(y m ) 1 is calculated all 
over the interference transmission range measured using 
the expression 

2hn {v m ) v m = m (1) 

The measurements were performed for plane-parallel 
plates of three thicknesses (from 0.8 to 1.2 mm) cut 
from silicon samples of all the isotope compositions using 
Bruker IFS-113v Fourier-transform vacuum spectrometer 
in 9500 . . . 110 cm -1 wavenumber range (1.05 ... 90 jim 



wavelength range) with 0.1 cm -1 resolution. The diam- 
eter of a radiation beam transmitted through the sam- 
ples did not exceed 3 mm. At least twenty data points 
were obtained for each interference maximum. The plates 
thickness was measured on DG-30 holographic submi- 
crometer with the accuracy better than 0.05 um. During 
transmission spectra and thickness measurements sample 
temperature was hold to be 21.2 ± 0.05 °C. 

Fig. [5] shows the refractive index spectral dependences 
over the entire measurement range for all the single crys- 
tal samples studied. As well the most frequently cited 
literature data for silicon single crystals with a natu- 
ral isotopic composition jl7H2ll | are shown. As evident 
from Fig. [2] the scatter of the refractive index data for 
silicon with a natural isotopic composition exceeds con- 
siderably the difference in the refractive index values of 
mono-isotopic silicon single crystals. 

With silicon isotope atomic mass increasing, the re- 
fractive index spectral dependences measured by us far 
from electron and phonon transition edges shifts down- 
wards. As an example, the refractive index values of the 
silicon single crystals are listed in Table [TTI for 1.5, 2 and 
5 jim wavelengths. 

It is known that two-phonon lattice absorption bands 
at 16.4 um (610 cm -1 ) with the maximum inten- 
sity of 9 cm -1 are observable in the 6. 7... 20 um 
(1500... 500 cm -1 ) spectral range in the spectra of 
1 . . . 2 mm-thick silicon samples. According to [22j, those 
should give rise to well-distinguished features in the re- 
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Figure 2: Spectral dependence of the refractive index: 
our measurments for silicon samples (a) 28 Si, (b) nat Si, 
(c) 29 Si, (d) 30 Si; the literature data for silicon crystals 
with natural isotopic composition (i) [17j . (h) [18| . 

(g) m, (e) m, ( f ) m 



Table II: Refractive index values of silicon single 
crystals at the wavelengths 1.5, 2, and 5 um 



Single crystal 


1.5 um 


2 um 


5 um 


28 Si 


3.48207 


3.45260 


3.42194 


natgj 


3.48191 


3.45243 


3.42177 


29 Si 


3.48171 


3.45224 


3.42161 


30 Si 


3.48112 


3.45168 


3.42107 



fractive index dispersion. The transmission spectra mea- 
sured with 2 cm -1 resolution for all 1.2 mm-thick sili- 
con crystals are shown in Fig. [3] (a). Fig. [3] (b) presents 
the variation of the crystals refractive index measured by 
us in this spectral region. The absorption band shift to 
longer wavelengths with silicon isotope atomic mass in- 
creasing and the anomalous behavior of refractive index 
dispersion well-conforming to the classic Lorentz theory 
are clearly seen. As far as we know, we managed to ob- 
serve such a refractive index spectral behavior in the re- 
gion of so low absorptionfless than 10 cm -1 ) for the first 
time. For example, in [20j (line e in Fig. [2]) the interfer- 
ence refractometry method was used as well to measure 
the refractive index of n-type silicon (phosphorus-doped 
with 3 ... 4 x cm resistivity) in 2.43 ... 25 um wave- 
length range at 26 °C. However, no peculiarities in the 
two-phonon absorption region were revealed. In the first 
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Figure 3: (a) Transmission spectral dependence of for 
1.2-mm thick 28 Si, 29 Si, 30 Si, and nat Si samples in the 

IR region of phonon absorption peak, 
(b) Refractive index spectral dependence of the same 
single crystals. 



approximation, the uncertainty in the refractive index 
values was attributed by the authors of 20] to the un- 
certainty in sample thickness (1 part per 10 4 ), which, in 
our opinion, is quite enough to observe the changes in 
the spectral dependence of silicon refractive index shown 
in Fig! (b). 

The refractive index values calculated from the trans- 
mission spectra are conventionally approximated by ana- 
lytical formulas allowing one to preserve the experimen- 
tal accuracy of refractive index determination. The Sell- 
meier or Herzberger formulas are often used to approx- 
imate the refractive index spectral dependence. How- 
ever, our studies show that the application of generalized 
Cauchy dispersion function 23] allows us to reduce the 
root-mean-square deviation of the approximation from 
the calculated refractive index values by several times. 
Therefore we used it to fit the experimental data: 
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Figure 4: Spectral dependence of the difference between 



the calculated, 



and approximated, n a 



refractive index for Si single crystal 



Here v is the wavenumber, v [cm *]= 10 4 x A 1 [pa]. 

Due to the presence of phonon absorption lines, the 
range of refractive index approximation was subdivided 
into four parts. This allowed us to reduce considerably 
the root-mean-square deviation between the measured 
and approximated refractive index values and to bring 
it to the level not exceeding 1 x 1CP 4 . The coefficients 
of Cauchy approximation ([2]) for refractive index of all 
silicon single crystals corresponding to the best fit to the 
experimental results, are given in Tables IIIIHVI1 The 
values of mean-square deviation for all approximation 



ranges are listed at the foot of the tables. As an ex- 
ample, the spectral dependence of the difference between 
the calculated refractive index for 30 Si single crystal and 
its approximated values is shown in Fig. 0] It is seen 
that the deviation of the approximation from the calcu- 
lated refractive index values does not exceed 5 x I0" 5 in 
a sufficiently wide range. However, with the wavelength 
exceeding 16 urn the accuracy of peak position determi- 
nation decreases owing to noticeable growth of noise in 
the transmission spectra. Accordingly, the accuracy of 
the refractive index calculation is reduced. 

Comparing our results to the data published, we found 
that the refractive index values closest to ours for natural 
silicon are given in [2ll | . In that work measurements were 
carried out by the prism method with I x I CP 4 relative ac- 
curacy at several wavelengths in the 1.1 . . .5.5 urn range 
at temperature from 20 to 300 K. The results obtained 
were approximated by the Sellmeier formula. 

B. On the origin of the dispersion curve shift 

The observed shift of the dispersion curve is easily un- 
derstood using classic models of a lattice of oscillators 
with dipolar coupling [3, Neglecting the spatial 

dispersion the dielectric constant of a cubic lattice is ex- 
pressed in this model by the equation 



1 



OLiEg {Mi 



fan (Mi) 
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where Mi is the atomic mass of J Si isotope, E g (Mi) is 
the band gap (fundamental absorption short-wavelength 
edge energy), f2 (Mi) is the phonon frequency (fundamen- 
tal absorption long- wavelength edge energy), on and Pi 
are the model parameters representing two types of oscil- 
lators corresponding to the electron and phonon funda- 
mental absorption edges, respectively. With no absorp- 
tion, the refractive index equals to n (oj) — e 1 ! 2 (oj). 

The fundamental absorption edge energies depend on 
isotopic composition and for an isotopically pure crystal 
can be expressed as (26l - [28j : 



E g (T, Mi) = E B -a B 



(M n 



n (Mi) = n 



(Mn 



V Mi 

1/2 



1/2 



V Mi 



exp(ft (M l )/T)-1_ 



where M nat is the average atomic mass of natural sil- 
icon, {exp (O (Mi) /T) — 1} is phonon Bose distribu- 
tion function, £1q is the model parameter. According to 
H|, E B w 28713.3 cm- 1 , a B « 967.9 cm" 1 . 

The oscillator parameters, a$ and fy, depend on iso- 
topic composition as well. First, both parameters depend 
on the lattice constant as a^~ 3 and ao is known to change 
with isotopic composition (see e.g. [29j]). However, this 



dependence is negligible in the qualitative model under 
consideration. Second, in the case of silicon the phonon- 
related parameter, Pi, describes the IR absorption due to 
two-phonon transitions and hence depends on the atomic 
mass as fa ~ Mr 2 [30]. Therefore 

on = a Q 

Pi = PoM~ 2 



Table III: Parameters of the approximation of 28 Si single crystal refractive index 





Spectral range 


Parameter 


9400... 745 cm" 1 
(1.06 . . . 13.42 pm) 


745 .. . 625 cm" 1 
(13.42 . . . 16.00 pm) 


625... 583 cm" 1 
(16.00 . . . 17.15 pm) 


583... 117 cm" 1 
(17.15 . . . 85.47 pm) 


A 6 


4.Uo4oo X 1U 


u.u 


u.u 


A A 71 90 v 1 ft - 
— 4.4f izy X 1U 


A A 


1.05838 x 10" 18 


0.0 


8.07309 x 10" 11 


3.89667 x 10" 14 


A 2 


1.42643 x 10" 9 


-1.75666 x 10" 6 


-1.18530 x 10" 4 


-1.00590 x 10" 8 


C 


3.41621 


6.74524 


68.44582 


3.41782 


B 2 


-34.48103 


-2.35655 x 10 6 


-1.58031 x 10 7 


-44.46299 


Bi 


0.0 


7.39436 x 10 11 


1.43557 x 10 12 


0.0 


B 6 


0.0 


-8.68186 x 10 16 


0.0 


0.0 


R 


2.0 x 10" 5 


3.0 x 10" 5 


9.0 x 10" 5 


1.1 x 10" 4 


Table IV: Parameters of the approximation of 29 Si single crystal refractive index 




Spectral range 


Parameter 


9400 ... 730 cm -1 
(1.06... 13.70 pm) 


730 .. . 617 cm" 1 
(13.70... 16.21 pm) 


617... 573 cm" 1 
(16.21 . . . 17.45 pm) 


573. . . 114 cm" 1 
(17.45 . . . 85.72 pm) 


A 6 


4.16405 x 10" 27 


0.0 


0.0 


0.0 


A A 


1.03238 x 10" 18 


0.0 


1.45163 x 10" 10 


4.61074 x 10" 15 


A 2 


1.42694 x 10" 9 


-1.35827 x 10" 6 


-2.03966 x 10" 4 


-5.60442 x 10" 10 


C 


3.41588 


5.89678 


110.64230 


3.41637 


B 2 


-29.06141 


-1.69322 x 10 6 


-2.49974 x 10 7 


7.04348 


Bi 


0.0 


5.12630 x 10 11 


2.18071 x 10 12 


-7.25535 x 10 s 


B 6 


0.0 


-5.81255 x 10 16 


0.0 


0.0 


R 


2.0 x 10" 8 


3.0 x 10" 5 


1.0 x 10" 4 


1.0 x 10" 4 


Table V: Parameters of the approximation of 30 Si single crystal refractive index 




Spectral range 


Parameter 


9310 ... 720 cm -1 
(1.07... 13.89 pm) 


720 .. . 605 cm" 1 
(13.89... 16.53 pm) 


605... 565 cm" 1 
(16.53 . . . 17.70 pm) 


565... 112 cm" 1 
(17.70 . . . 89.29 pm) 


A w 


0.0 


1.58348 x 10" 28 


0.0 


0.0 


As 


0.0 


-3.70584 x 10" 22 


0.0 


0.0 


A 6 


3.87178 x 10" 27 


3.45410 x 10" 16 


0.0 


-5.25168 x 10" 21 


A 4 


1.06782 x 10" 18 


-1.60314 x 10" 10 


1.37221 x 10" 10 


1.67675 x 10" 14 


A 2 


1.42449 x 10" 9 


3.70631 x 10" 5 


-1.87828 x 10" 4 


-6.42102 x 10" 9 


C 


3.41535 


0.0 


99.59448 


3.41676 


B 2 


-50.69159 


0.0 


-2.18377 x 10 7 


-43.76842 


Bi 


0.0 


0.0 


1.85529 x 10 12 


0.0 


R 


2.0 x 10" 8 


3.0 x 10" 5 


5.0 x 10" 5 


1.0 x 10" 4 


Table VI: Parameters of the approximation of nat Si single crystal refractive index 




Spectral range 


Parameter 


9480 ... 745 cm -1 
(1.06... 13.42 pm) 


745 .. . 625 cm" 1 
(13.42 . . . 16.00 pm) 


625... 583 cm" 1 
(16.00... 17.15 pm) 


583. . . 125 cm" 1 
(17.15... 80.00 pm) 


A 6 


3.85272 x 10" 27 


0.0 


0.0 


0.0 


A A 


1.08707 x 10" 18 


0.0 


1.06100 x 10" 1() 


1.10503 x 10" 14 


A 2 


1.42510 x 10" 9 


-1.38069 x 10" 6 


-1.54334 x 10" 4 


-4.28163 x 10" 9 


C 


3.41618 


6.01658 


87.37999 


3.41725 


B 2 


-96.96621 


-1.82986 x 10 6 


-2.02501 x 10 7 


-51.15422 


Bi 


0.0 


5.71062 x 10 11 


1.82698 x 10 12 


-1.93263 x 10 s 


Be 


0.0 


-6.67343 x 10 16 


0.0 


0.0 


R 


2.0 x 10" 8 


2.0 x 10" 5 


6.0 x 10" 5 


1.3 x 10" 4 
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Figure 5: Spectral dependence of the left side of 
equation [3] for optimal and differing by ±1 values of 
mi o. 



where a$ — const, /3q = const. 

Thus, with Si atomic mass increasing, the electronic 
absorption edge shifts to shorter wavelengths, and the 
phonon absorption edge shifts to longer wavelengths. As 
a result, the dispersion curve of the refractive index in 
the transmission window should shift downwards with- 
out considerable changes in shape, as observed in our 
experiments. 



VI. ACCURACY ANALYSIS 

To measure the refractive index precisely by interfer- 
ence refractometry, the orders, m, of interference max- 
ima should be determined with absolute accuracy. For 
this purpose we used the relation following from expres- 
sion (JTJ) written for two samples of different thickness, 
hi and hi, with interference orders mi (v) and m 2 (v), 
respectively: 

mi M hi 

— — = const (3) 

m 2 (v) h 2 

If sample thickness is known with absolute accuracy, 
the constant in the right side of this equation turns to 
zero for a certain relation between mi (y) and m,2 (v). 
However, usually the sample thickness can be measured 
with some finite accuracy resulting in non-zero value of 
this constant [l6| . Moreover, if the transmission spectra 
of the samples of different thickness are measured at dif- 
ferent temperatures, spectral dependence of the left side 



of equation (j3]) will be affected by the temperature de- 
pendence of the refractive index of a material. However, 
if one makes a wrong choice of the maxima orders at least 
by 1 , the spectral dependence of the left side of equation 
([3]) will explicitly differ from constant as evident from 
Fig- for two nat Si samples. 

Minimization of the right side of dependence ([3]) allows 
us to determine uniquely the maxima orders for each pair 
of samples. Additionally, the procedure of interference 
orders determination becomes simpler if the formula ([T]) 
is considered to be valid not only for the maxima but as 
well for all other points of the spectral transmission curve 
for which m (y) values are not integer. One should use 
a linear approximation for mi. 2 frequency dependence in 
the interval between two nearest maxima (according to 
our calculations such an approximation is valid with the 
accuracy about 1CT 6 for 1 ... 3 cm distance between 
the maxima) and increase the corresponding order by 1 
passing over the next maximum. In this case the ap- 
plication of the ([T]) criterion becomes a simple search of 
orders corresponding to a given frequency in the trans- 
mission spectra of sample pairs, since the set of frequency 
values is the same for all spectra. 

In the method of interference refractometry the error of 
the refractive index calculation depends on determination 
errors of all the variables occurring in the formula ((T|). 
If Fourier-transform spectrometer with the wavenumber 
precision better than 0.01 cm -1 is used, this error is 
mostly contributed by the sample thickness measurement 
error, Ah. The corresponding error in silicon refractive 
index is equal to An = —nAh/h. If Ah — 0.05 um and 
values of m are determined precisely for all interference 
maxima, one gets Art k 2 x 10~ 4 on average for the only 
~ 1 mm-thick sample used. 

Using of two and more samples allows us to reduce 
this error due to the correlation of the refractive index 
values determined for all samples over the entire mea- 
surement range. Assuming the refractive index values 
to be identical for all the samples within the entire mea- 
surement range, we can analyze the influence of accuracy 
of the used thickness values on the final result of mea- 
surements and "correct" these thickness values. Analysis 
of the discrepancy between the "corrected" and initially 
used thicknesses provides an idea of how well the sam- 
ples are prepared and measurements are made. Correc- 
tion of thickness values (hi, i = 1,2,3 in the case of 
three samples) is performed by orthogonal projecting in 
the (hi, /12, /13) space the point with coordinates equal 
to these thickness values on a straight line representing 
solution of the following set of equations for all pairs of 
samples 



hi mi hi mi 

hi m 2 h 3 m 3 



ho mo 

0, T 1 - — 

h 3 m 3 







The refractive index values scatter, Am, for sample of 
thickness hi is obtained immediately from the relation 



An,: 



Ah, 



Wavelength, urn 

1 10 100 

1 1 1 1 1 — i — i — r~| 1 1 1 1 1 — r~ri — r~ 




10000 1000 100 

Wavenumber, cm" 1 

Figure 6: Spectral dependence of the difference Ani23 
averaged between refractive index for 28 Si samples of 
different thicknesses. 



where n$ is the "averaged" refractive index value, Ahi is 
the difference between the "corrected" and initial thick- 
ness. 

Table IVIII lists sample thickness measured on submi- 
crometer at 21.2 °C and the thickness values corrected 
using the interference spectra in the middle (" MIR" , from 
2 to 25 urn) and near ("NIR" , from 1.06 to 2 pm) infrared 
ranges, as well as the corresponding difference between 
the calculated and measured thicknesses. As is seen from 
Table IVIII the thickness correction value turns out to be 
extremely small for most of the samples not exceeding 
the error of direct measurement of the thickness. We es- 
timate the average error of the refractive index determi- 
nation resulting from the error of the samples thickness 
determination to be 1.3 x 10~ 4 . The contribution of the 
remaining factors influencing the samples transmission 
spectra measurements and, accordingly, our calculation 
of the refractive index performed similarly to paper [3l| , 
is found to be 10 times less and therefore is negligible. 

Unlike the errors of thickness measurement and inter- 
ference maxima orders determination, the error of the 
maxima positions determination is the only one exhibit- 
ing distinct frequency dependence: the latter error grows 
with the frequency decreasing (wavelength increasing) or 
in the absorption band vicinity and near the transparency 
range edges if any in the transmission spectra. To deter- 
mine the maxima positions more exactly we restricted 
the range of the sample thicknesses used and spectrome- 
ter resolution to get not less than 20 experimental points 
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Table VII: Measured and corrected thicknesses of Si 
samples (in pm; i — 1, 2, 3) 





hi (measured) 


hf 11 (MIR) 


hf" (NIR) 


A/lMIR 




28q- 


815.87 

i nnQ ok 

1179.09 


815.98 
1179.03 


815.99 

i nnQ oq 
1179.02 


0.11 

n no 
— U.Uz 

-0.06 


0.12 

n no 
— U.Uz 

-0.07 


29 Si 


815.95 
1003.23 
1179.24 


815.98 
1003.21 
1179.23 


815.98 
1003.23 
1179.22 


0.03 
-0.02 
-0.01 


0.03 
0.0 
-0.02 


30 Si 


816.09 
1003.34 
1179.54 


816.12 
1003.41 
1179.46 


816.13 
1003.38 
1179.48 


0.03 
0.07 
-0.08 


0.04 
0.04 
-0.06 


natgj 


815.87 
1003.27 
1179.15 


815.88 
1003.28 
1179.14 


815.88 
1003.28 
1179.14 


0.01 
0.01 
-0.01 


0.01 
0.01 
-0.01 



between adjacent maxima in the corresponding transmis- 
sion spectra. To estimate the influence of error of the 
maxima positions determination on the refractive index 
calculation, we determined the averaged difference be- 
tween the refractive index values calculated for samples 
with different thickness. As an example, such a difference 
is presented in Fig. [5] for 28 Si samples. 



VII. SUMMARY 

In this work we have prepared and characterized high- 
purity single crystals of silicon stable isotopes, 28 Si, 29 Si, 
and 30 Si with enrichment above 99.9 at.%, and of sil- 
icon with natural isotopic composition, nat Si. For the 
first time the spectral dependences of the refractive index 
of all single crystals are measured in wide transparency 
range from 1.06 to about 80 pm including the phonon 
absorption region. 

The results of our study prove the refractive index val- 
ues of mono-isotopic silicon single crystals to differ con- 
siderably. This may provide a possibility of manufacture 
of multilayer wave-guiding and quantum-sized structures 
operating in the IR spectrum range from different silicon 
isotopes to solve various problems of electronics and fiber 
and integrated optics. 
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